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We present a scheme that transform 1 qubit to A4 identical copies with optimal fidelity via free 
dynamical evolution of spin star networks. We show that the Heisenberg XXZ coupling can fulfill 
the challenge. The initial state of the copying machine and the parameters of the spin Hamiltonian 
are discussed in detail. Furthermore we have proposed a feasible method to prepare the initial state 
of the copying machine. 
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One of the most fundamental difference between clas- 
sical and quantum information is the no-cloning theorem 
[1,2]. It states that accurate cloning of any arbitrary 
quantum state is impossible. Nevertheless it doesn't for- 
bid one to clone quantum states approximately. In the 
early work of Buzek and Hillery an optimal 1—^2 
universal quantum cloning machine (UQCM) was pro- 
posed where the copying process is input-state inde- 
pendent. And quantum cloning machine for equatorial 
qubits which are so called phase-covariant cloning (PCC) 
machine was proposed by Brufi et al. [1] . Optimal fidelity 
and optimal quantum cloning transformation of general 
TV to M (M > N) are presented in [5-8]. It was also 
shown a UQCM can be realized by a network consisting 
of quantum gates 0. 

Several approaches have been made to realize the 
unitary transformations leading to the cloning process 
experimentally [10-13]. However, most of these schemes 
are based on quantum logic gates and post-selection 
methods, which need time modulations. Recently, quan- 
tum computation via spin networks based on Heisenberg 
couplings was presented [14-27]. One achieve is that with 
Heisenberg chains, high fidelity quantum state transfer 
can be achieved [15-24]. The most attracting feature of 
this approach is that it needn't time modulation for the 
qubits couplings. Once the initial states and the evolu- 
tional Hamiltonian is determined, the system can faith- 
fully implement designated computation task through 
free dynamical evolution. The whole computational evo- 
lution does not involve any external controlling, which 
provides relatively longer decoherence time for the sys- 
tem. Schemes for PCC via spin networks was proposed 
in the work of De Chiara et al. [2^, [2^ . Chen et al. [27| 
further improved the 1 ^ M PCC case to an optimal 
level. However the optimal UQCM via a spin network is 
still a challenge. 

In this paper, we show that by properly introducing 
the ancilla qubies, designing the spin exchange interac- 
tions, and choosing the initial state of the cloning ma- 
chine, optimal 1 ^ M UQCM can be realized via the 
free evolution of a spin star network Hamiltonian. More- 
over a scheme on preparing the initial state of the cloning 
machine have been proposed. 



The spin network involved in our scheme forms a star 
configuration (See Fig. 1(1)). The central qubit (input 
state) is labeled /, the M target qubits labeled T, and 
the M — 2 ancillas labeled A. We start with the conven- 
tional Heisenberg XXZ coupling Hamiltonian without an 
externally applied magnetic field. 
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where aj'^'^, cr^:^'^, a^^'^ are Pauli matrices of the input 
particle, the target qubit, and the ancilla qubits respec- 
tively (we introduce M — 2 ancilla qubits), Ji and J2 are 
the exchange spin coupling coefficients between the in- 
put qubit with the target qubits and the ancilla qubits 
respectively, Ai and A2 are the anisotropy parameters 
(when A = 0, the Hamiltonian reduces to XX model 
while A = 1 it corresponds to Heisenberg model). 

Following with Gisin and Massar 0] we suppose the 
unitary transformation for optimal 1 —^ M cloning take 
the form: 
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where Ui^m — e^'^*" (tQ is the evolution time) denotes 
the free evolution of the spin system, \R) denotes the ini- 
tial state of the copying machine and AI blank copies. 
\S{M,i))T is the normalized symmetry state of the M 
target qubits with i spins up. \Ri) are orthogonal nor- 
malized sates of the ancilla qubies (here include the input 
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FIG. 1: (1) Spin star network for 1 ^ M UQCM. The left 
side spins form the ancilla, and spins on the right side form 
the target particles. (2) Spin network employed for generating 
initial states. Each ancilla spin interacts with all the target 
spins. 



qubit). We choose the initial state \R) as follows: 



\R) = cj2 vm - 

|a,) = \SiM,i))T<i^\S{M-2,M-l-i))A, 
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where C 



is the normalization factor, 
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\S{M — 2, Af — 1 — i))A is the normalized symmetry state 
of the M — 2 ancilla qubits. Noticing \R) is invariant 
under the spin flipping operation, we will show later by 
spin flipping both sides of Eq.([21), Eq.Q is automatically 
satisfied. 

First we discuss the conditions to satisfy Eq.®. In- 
stead of studying all the states in the Hilbert space of 
the Hailtonian ([1]), we would rather to introduce a two 
dimensional subspace Hat (we use ^pab to note states in 
this subspace, V'ab to note states orthogonal to this sub- 
space) , which is spanned by two basic normalized orthog- 
onal states I a), and \b) 



\a) = I T)/® , 
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\b,) = \S{M,j + 1))t®\S{M-2,M-1-j))a ■ 

Notice I a) is our initial state for Eq.([2]), and we will show 
that some linear combination of these two states has the 
same form of the righthand side of Eq. ^ . We find if the 
parameters of the spin Hamiltonian ([1]) obey the follow- 
ing relations 
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the subspace we choosing is closed, i.e., (V'a&I^IV'ofc) = 0- 
Then it is convenient for us to calculate the free evolution 
of the system in this two dimensional subspace. It is 
useful to rewrite the Hamiltonian ([T]) with the Ladder 
operators. Using the relations ([7]) the Hamiltonian take 



the form: 

S± = K±Zaf)/2, 4=^4: ji=T.'A, 

T A 

s| = CTo/2, = '^a^/2, J\^^a\/2, 
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H - j{s+{J^-J2) + s^{J+-J+) + 2\s^jm + J%)). 

With this representation of the spin Hamiltonian it is 
easy for us to calculate H act upon our bases. 

H\a) - V2^|6), 

H\b) = -J\\b)+V2J\a). 

Thus, the Hamiltonian is closed in the subspace Tiab, and 
we can write the matrix form of H in Tiab, 
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This is the key point of our scheme. Despite how large M 
is, such a two dimensional space always exist as long as 
the condition ([7]) is kept. Now, our problem reduced to a 
two dimensional quantum evolution in TLab^ the unitary 
transformation Ui_M{t) takes the form: 

Ui,Mit) - exp{-ii/t} 

-1 ^xTZps (^'^^ + 2 V2a,)) , 

where tr^ and are pauli matrices in Tiab- We choose 
the anisotropy parameter 
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As our initial state for Eq.® is |a), after having evolved 
for t. 



Ht)) = Ui,M{t)\a) 
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where ip — \pijt is the rescaled time parameter. When 
ip = 7r/2, i.e, 



t = tQ = 



(9) 



the state of the system take the following form. 
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where \S{M — 1, i))A<gii denotes the normahzed symmetry 
state in the direct product space of the input qubit and 
ancilla qubits. The state \a{to)) is exact the same form as 
the righthand side of Eq.©- The orthogonal normahzed 
states \Ri) take the form: 



|5(M-l,z))^«/. 



To go further, we introduce the spin flipping operator, 



M 
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This unitary operation flip aU the spins in our consid- 
eration. It is easy to see that the Heisenberg XXZ spin 
Hamihonian ([1]) is invariant under such operation, i.e., 
PHP-^ = H. The initial state for Eq.© is | i)/ ® \R) = 
P I a) = I Op), after having evolved for t 




|ap(t)) = f/i,MW|ap> - P Ui,M{t)\a) = P \a{t)). 



When the evolution time t ~ to 



|ap(io)) = P Hto)) = -ie'-^* I^(A^'*))t 

® \S{M ~l,M ^1-i))a<^i. 

It is exact the same form as the righthand side of Eq. ^ . 
The above calculation show that we can find such condi- 
tions ()4|7l8|9p satisfying Eq.(l2l[3|) simultaneously, i.e, the 
optimal cloning can be fulfilled under such conditions. 

One interesting thing is that through the beginning 
to the end of this free evolution the fidelity of a single 
copy to the input is independent of the input state (a 
universal cloning) . Suppose the input state is: | input) / = 
Q^l T)/ + /3| After having evolved for i, the state 
of the system take the form: \t) — oi\a{t)) + P\ap{t)). 
The reduced density matrix of a single copy at t can be 
calculated directly, 
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The fidelity of this copy is 
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F is only a function of the rescaled time ip {(p = VSJ^t). 
So the whole cloning process is input state independent. 
When t = the fidelity is 1/2, and when t — to the 
fidelity reaches it's optimal bound {2M + l)/3Af (see 
Fig.2). 

One shortcoming of quantum cloning based on logic 
gates is the circuit becomes more complicated as M in- 
creases. As a result, when M is large it may be difficult 



FIG. 2: The input state independent fidelity of a single copy 
as a function of the rescaled time ip {ip — VSJ^t), for M = 
2,3,4. When ip = nn the fidelity equals 1/2, and when p — 
(2n 4- l)7r/2 the fidelity reaches it's optimal bound {2AI + 
1)/3M. 



for one to go through the copying process before the state 
having been decoherenced. However, the evolution time 
of our scheme is = \/37i'/(12j/), which is independent 
of M. This is an advantage to fulfill large M cloning. 

The problem now is how to prepare the initial state 
For M — 2 (M = 3), |i?) is two (four) particle symmetry 
state. But for M > 3, \R) is not simply a symmetry state. 
Interestingly, we find that \R) is exactly the ground state 
of some spin Hamiltonian. And it is feasible for one to 
prepare it by just cooling the system. Such Hamiltonian 
is consisted of two parts 



H' = H'a + H[. 



(11) 



Hq is the part with Heisenberg XXZ coupling (A — 
— 1) between the target qubits and the ancilla qubits 
(Fig.l(2)), 



(12) 



where J' is the spin coupling coefficient, Jt and J a are 
total angular momentum operators of the target qubits 
and the ancilla qubits respectively. H[ is the part with 
Ising coupling between all the qubits, 
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where A is the coupling coefficient. These two parts are 
commute, [Hq, H[] — 0. We find \R) is an eigenvector of 
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Hq and H[ simultaneously 

To prove \R) is the ground state we solve the spectrum of 
H'. We introduce the unitary operator Qt to act on the 
target qubits (it is equivalence to introduce Qa acting on 
the ancilla qubits), 

M 

Qt = Qt^ ~ H '''Ti • 

1=1 

This unitary operation transforms Hq to the Heisenberg 
Hamiltonian and leaves H[ unchange, 

QtHqQj.^ = —J^'{J^J^ + + 2 J^) 

So the spectrum of H' is 

E' = -J'(j(j + 1)-Jt(jt + 1)-Ja(ja + 1)) 
+A(j-)2 - 



^.,,2 A(M - 1) 



where j (j^) is the total {z component) angular momen- 
tum quantum number of the transformed Hamiltonian. 



If we choose J' < 0, and A > 0, the nondegenerate 
ground state energy of H' is J'{M'^ - 4)/2 - A(M- l)/2 
(j = l,jT = M/2,jA = (M - 2)/2,r = 0), which is 
just the eigenvalue of \R). So far, we have proved \R) 
is the ground state of H' . Thus the initial state of the 
copying machine can be prepared by cooling the system. 
No measurement is involved in this implementation, and 
also we needn't any time modulation of the Hamiltonian. 



Through out this paper, optimal UQCM that produce 
M copies out of a single input via free evolution of spin 
star networks has been discussed. We have proved for 
arbitrary M the unitary evolution can be fulfilled in a 
two dimensional subspace. Using this character we find 
the analytical solutions for the optimal 1 — > Af univer- 
sal cloning process. Through this process the fidelity 
keeps input state independent, and it reaches the op- 
timal bound at t — \/3n/6J^, which is independent of 
M. Also we have studied the initial state of the coping 
machine in detail, and find it is exactly the ground state 
of some spin Hamiltonian (only quadratic terms are in- 
volved). Thus, the preparation of the initial state can 
be accomplished by cooling such systems. No measure- 
ment and time modulation is involved here. Therefore 
our result opens up a promising prospect towards robust 
optimal UQCM. Such a p rospect is relevant for several 
experimental systems j28l f29l|. 
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